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In 0,0 I. M. Suslov proposes a theoretical description for the interplay between disorder and 
superconductivity that, among other things, claims to predict situations where superconductivity 
is enhanced by disorder. In this comment we show that Suslov’s results do not make any sound 
predictions relating to this problem. Suslov also suggests that the percolation approach recently 
employed by the authors to compute the global critical temperature of inhomogeneous superconduc¬ 
tors induced by weak multifractality 0 is not satisfactory. We refute Suslov’s opinion using simple 
ideas from the Bardeen-Cooper-Schiaffer (BCS) theory of superconductivity. More importantly we 
stress that our theory makes quantitative analytical predictions for observables such as the global 
critical temperature and the spatial distribution of the order parameter. Our work can therefore be 
easily confirmed or disproved through either numerical or experimental study. 
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Suslov has recently claimed 0 to have developed a the¬ 
ory, for the interplay of superconductivity and disorder, 
that is capable of answering questions like: Can disorder 
enhance superconductivity? He also criticised [l| the per¬ 
colation approach proposed by the authors Q to address 
this problem in the case of metallic superconductors with 
inhomogeneity induced by weak multifractality. In the 
following we refute this criticism and show that Suslov’s 
theory is based on an ill defined critical temperature and 
the unjustified use of the mean-field approximation. Un¬ 
surprisingly, Suslov’s main prediction of a large enhance¬ 
ment of the critical temperature driven by disorder has 
never, and very likely will never, been found numerically 
or experimentally. 

Suslov studies the critical temperature of a supercon¬ 
ductor T c after adding a small concentration of point-like 
impurities in two limits: 

1. Strong coupling: corresponding to an impurity 
that act as a deep potential well of the size the 
lattice spacing, a, and is effectively isolated from 
the rest of the system. This region is completely 
unrelated to our recent investigation 0. 

2. Weak coupling: corresponding to a weak random 
potential where the superconductor is not strongly 
affected by the impurity. 

Let us first discuss the case of strong coupling. Suslov as¬ 
sumes that these impurities induce discrete levels around 
the Fermi energy and then uses mean-field techniques to 
compute the critical temperature. His main result is that 
T c ~ E 0 A where E 0 is of the order of Fermi energy and 
A is the electron-phonon coupling constant. However, 
it is well established that a mean-field approximation is 
completely unjustified in this ultra-small limit where the 
effective mean-level spacing <5, controlled by the lattice 
spacing a, is much larger than the bulk gap Ao- 

In this region, rightly termed catastrophic by Suslov, 


the spectroscopic gap must be computed using either the 
Richardson equation 0] or perturbative techniques 0] 
starting from the normal state. More importantly, the 
critical temperature T c is ill defined in this fluctuation 
dominated regime 0 where the spectroscopic gap at the 
Fermi energy is always non-zero at finite temperature 
0 0 • Even for Ao ~ 6 00 there is no substantial decay 
of the spectroscopic gap at the mean-field critical temper¬ 
ature. This is not surprising since for sizes comparable 
to the lattice spacing, superconductivity is only a weak 
perturbation 0. Suslov’s prediction of T c ~ 5000K, for 
values corresponding to typical metallic superconductors, 
is therefore completely, though not surprisingly, unreal¬ 
istic. 

Suslov goes on to argue, reasonably, that since this 
critical temperature only occurs in a few, well sepa¬ 
rated regions, of typical size given by the lattice spac¬ 
ing a, it cannot be the true global critical temperature 
of the sample. More specifically this T c is not the one 
at which global phase-coherence, and zero resistance, oc¬ 
curs. Suslov claims that it is possible to tune the number 
of resonant impurities, controlled by T c /Eq , by simply 
increasing the density of impurities, so that, 11 eventually 
their Josephson interaction is strong enough for stabiliza¬ 
tion of the mean-field solution at practically the same T c 
value". We first note this heuristic picture of emergent 
granularity for establishing phase-coherence in a disor¬ 
dered superconductor was proposed a long time ago 00 
and is still actively pursued [10]. However, the way Suslov 
implements it is open to criticism. First, there is broad 
evidence that this appealing picture is correct starting 
from the insulator side and gradually reducing disorder 
so that phase coherence is eventually achieved. Surpris¬ 
ingly Suslov applies it in the opposite direction since, 
according to him, phase coherence occurs by gradually 
increasing disorder. It is unclear to us how any genuine 
realisation of disorder can have this property. Disorder 
makes the system more localized and therefore always 
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increases the magnitude of phase fluctuations in the sys¬ 
tem. We stress that as disorder decreases from the deep 
insulator limit superconducting spots start to appear but 
are related to relatively flat regions of the disorder not 
to any resonant scattering 

Even assuming Suslov’s picture and its prediction of 
T c to be accurate, which we have shown above is not the 
case, in a weakly coupled superconductor A <C 1 so T c /E 0 
is always small and the global critical temperature must 
necessarily be zero. Moreover we have serious reserva¬ 
tions as to whether it is possible to tune this number in 
any realistic disordered material. In summary, Suslov’s 
statement that, after including phase fluctuations, the 
global critical temperature is still not very different from 
T c ~ 5000K is absolutely unrealistic. 

Towards the end of [l[ Suslov argues, and we agree, 
that Coulomb interactions will eventually reduce the 
global critical temperature to a value much smaller than 
T c ~ E 0 A <C 5000K. Unfortunately he does not give any 
quantitative or semi-quantitative estimate of the suppres¬ 
sion of the critical temperature due to this effect. The 
suppression of T c will strongly depend on the typical size 
of the superconducting drops and the coupling strength 
between neighbouring grains. There is recent experimen¬ 
tal evidence that even for almost isolated superconduct¬ 
ing nano-grains of Sn and Pb El of size ~ 7nm there 
is no Coulomb blockade or even substantial reduction of 
the superconducting gap. It is therefore unlikely that 
Coulomb interactions can reduce the critical tempera¬ 
ture by the several orders of magnitude Suslov’s the¬ 
ory implies, at least for metallic superconductors. This 
is another indication that Suslov’s critical temperature, 
T c ~ E 0 A, is ill defined and not directly related to a 
typical energy of any realistic superconductor. 

In summary the enhancement of superconductivity 
predicted by Suslov in the catastrophic regime is un¬ 
physical. Not surprisingly it has never been observed 
and, from our point of view, it will never be observed 
numerically or experimentally. We must also emphasise 
that this model is completely unrelated to the weakly 
coupled inhomogeneous superconductor induced by 
weak multifractality that we have recently investigated 

i- 


We now turn to the weakly coupled disordered region 
of Ref. p], [2j . According to Suslov the critical temperature 
is given by, 


ST C _ 1 f d UQUi(r) + izi(r) 2 

T c0 ~ XL d J CL r vl 


( 1 ) 


where v\ ( r ) is the deviation of the local density of states 
at the Fermi energy, z/p(r), from its unperturbed value 
i/o, A = gv o is the dimensionless coupling constant, T c o 
is the transition temperature in the absence of disorder 
and L d is “a volume for one impurity. The length, L , 
is poorly defined as it is unclear whether it is related to 
the system size. Similarly v\ (r) is not properly defined 
as a function of a small parameter. On physical grounds 


we speculate that its strength is controlled by the dimen¬ 
sionless conductance, go, so that the small parameter in 
the expansion is 1/goX. Irrespective of this, the simple 
dependence on A indicates that this is a first order purely 
perturbative result that describes the leading correction 
to the Anderson theorem. Suslov uses this expression as 
the starting point for his criticism of our recent paper 
Q, in which we compute analytically the global critical 
temperature of a disordered superconductor induced by 
weak- mult ifr act lity. 

Before we address, and refute, his criticism we sum¬ 
marise our main results. In Ref. [ 3 } we compute the su¬ 
perconducting order parameter analytically, including its 
energy dependence and statistical distribution in space. 
The spatial distribution of the order parameter is found 
to be always log-normal. By applying percolation tech¬ 
niques we determine the global critical temperature ne¬ 
glecting phase fluctuations. We find that it is enhanced 
with respect to the clean limit only for very weakly cou¬ 
pled superconductors. This enhancement seems robust 
as it persists even after the percolation threshold is in¬ 
creased in order to crudely estimate the effect of phase- 
fluctuations. It is important to note that in a disordered 
superconductor there are different critical temperatures: 
the global critical temperature corresponding to the max¬ 
imum at which a supercurrent can flow, the local criti¬ 
cal temperature at which the order parameter at a given 
point vanishes and the critical temperature related to the 
vanishing of the gap at the Fermi energy. For instance 
the first is the natural outcome of a transport measure¬ 
ment and the second is the natural outcome of a scanning 
tunnelling microscope experiment El- One of the con¬ 
clusions of our recent paper Q is that a sharp division 
between these distinct, but related, critical temperature 
is a fundamental requirement if one is to make substan¬ 
tial progress on this problem. 

The main criticism of Ref. 1] is directed at the per¬ 
colation approach we use to compute the global critical 
temperature. According to Suslov, since we are work¬ 
ing in the weakly disordered regime of Eq. ©> the lo¬ 
cal critical temperatures T c (r), defined as the tempera¬ 
ture at which the amplitude of the order parameter at 
r vanishes, depends only weakly on r. Consequently the 
value of T c resulting from percolation is similar to that 
obtained by averaging over T c (r) and therefore a per¬ 
colation approach is not necessary. This is not correct. 
Below we argue that Suslov’s opinion is based on some 
misunderstanding of basic features of our model and of 
Bardeen-Cooper-Schrieffer (BCS) theory as well. 

It is well known [uj that multifractality, a key ingre¬ 
dient in our results, induces strong correlations [12|, |l3( 
among eigenstates with energies within the Debye win¬ 
dow, cd- These enhance the spatial inhomogeneities 
of the order parameter even in the limit of weak- 
multifractality; corresponding to a large dimensionless 
conductance. Moreover, in BCS theory, the order pa¬ 
rameter depends exponentially on the coupling constant, 
A. As a consequence, small corrections to the coupling 
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constant, due to disorder or other perturbations, induce 
an exponential effect on the order parameter provided 
that the calculation is carried out self-consistently, typ¬ 
ically involving a re-summation of perturbative correc¬ 
tions. This is why even weak perturbations can induce 
strong modifications of the gap and critical temperature. 
We must also stress that there is substantial experimental 
0 , numerical 0, 0] and analytical [HI, [l7| evidence 
that supports this statement. 

The 1/A dependence of Eq. JTJ reveals a simple first 
order perturbative result that does not take into account 
the non-perturbative effect of the correlations between 
multifractal eigenstates. By contrast the combination of 
multifractality and a fully selfconsistent analytical cal¬ 
culation leads to a broad log-normal distribution for the 
spatial distribution of the superconducting gap, the as¬ 
sociated local critical temperature, as well as to a much 
more involved dependence of the global critical temper¬ 
ature on the coupling constant A. It is evident, see Fig 
4 of Ref. !3j, that such a broad distribution leads to im¬ 
portant variations of T c (r) at different points in space. 
Moreover we are interested only in weakly coupled su¬ 
perconductors in the limit of weak perturbations where a 
mean-field approach is still acceptable and the necessary 
truncation of the Hilbert space, by expressing the ma¬ 
trix elements as a function of the one-body eigenstates, 
is still a good approximation 00 . This fully justi¬ 
fies the computation of an upper bound for the global 
critical temperature by a percolation analysis of T c (r). 
Indeed our results, see Fig. 4 of Ref. 3, indicate the 
global critical temperature is vastly different, and more 
importantly lower, than any averaged T c (r) equivalent 
to T c in Suslov’s paper, even before any reduction due to 
Coulomb interactions. 

We also believe part of Suslov’s criticism of Ref. [3} is 
simply based on a poor understanding of our percolation 
approach. According to Suslov our percolation analysis 
is intended to describe fluctuations. This is not the 
case. In a first approximation our percolation analysis 
neglects completely phase fluctuations. It only gives 
the maximum temperature for which a finite amplitude 
of the superconducting gap permeates the sample. We 
then go on to provide a rough estimation of the effect of 
phase fluctuations by increasing the percolation thresh¬ 
old so that, even if superconducting regions permeate 
the sample, there is no long-range order. This is the 
correct approach when starting from a purely metallic 
state which becomes increasingly inhomogeneous with 
increasin g di sorder. Indeed there is recent experimental 
evidence 0 that the spatial distribution of the order 
parameter in a metallic quasi two-dimensional disor¬ 
dered superconductor is similar to the one we derived 
analytically. 


Finally we comment on other statements in Suslov’s 
papers that are not correct. Suslov claims that our re¬ 
sults are only valid in small and strictly two dimensional 
systems. However it is well established in the theory of 
disordered systems that weak multifractality also occurs 
in quasi two dimensional systems (thin films) and two 
dimensional systems with spin-orbit interactions. All of 
these systems are of direct experimental relevance. 

Suslov also states, referring to our calculation of the 
critical temperature, that ” there is no need to use the 
approximate equation (6) [equation ( 6 ) in [4j] when the 
accurate result [restated here as equation[T]] is available”. 
We emphasise that precisely the opposite is true. Suslov’s 
formalism is strictly perturbative with no sharp distinc¬ 
tion between the different critical temperatures that de¬ 
fine a disordered superconductor. Self-consistency is not 
reinforced and consequently the dependence on A is triv¬ 
ial. Eq. [1] could be a good starting point to re-derive our 
results. However this will require much more effort than 
Ref. [l] suggests. 

Suslov claims towards the end of [lj that we do not 
discuss phase-fluctuations and Coulomb interactions. We 
must point out that section VII from Ref. @| does already 
discuss this in some detail from both a qualitative and 
semi-quantitative viewpoint. Suslov also states that we 
claim Eq is always much larger than cd- We are inter¬ 
ested in weakly coupled metallic superconductors so in 
our case this limit is always applicable. However, we 
make no such claim regarding other materials. 

An important part of Ref. [l| is actually devoted to 
criticising other recent publications on superconductivity 
with multifractal eigenstates 0,01 that are not directly 
related to our work. For example, the critical tempera¬ 
ture at which the gap at the Fermi energy vanishes is, 
in these publications, T c ~ EqX 1 ^ with 7 related to the 
multifractal dimensions. This is an expression that, we 
are at pains to stress, is not recovered in our formalism 
in the limit of weak-multifractility. 

To summarise, Suslov’s claim to have explored conclu¬ 
sively the possibility for enhancement of superconductiv¬ 
ity due to disorder are clearly unfounded. His theory 
is based on ill defined concepts, the results are rarely 
testable and lead, in general, to unphysical predictions. 
Similarly, the criticism of the percolation formalism em¬ 
ployed in our paper Q is not warranted. Percolation is 
a powerful tool for an quantitative study of the problem 
of inhomogeneous superconductivity induced by weak- 
multifractality. 
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